import pandas as pd
import geopandas as gpd
from shapely.geometry import Point
import re
from tqdm.auto import tqdm
import contextily as cx
cx.set_cache_dir("tile_cache")
df = pd.read_excel("Pollution analyses Mapua.xlsx")
df.rename(columns={"Aluminuim": "Aluminium"})
df
| Sediment | Location | Total DDT Isomers | Dieldrin | Aluminuim | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.0075 | 0.001 | 16518.912655 | 1052.643599 | 42.463242 | 49.985737 | 273.474991 | 8.104016 | 52.778163 | 8.868234 | 35.506181 | 5.43789 | 0.991537 | 0.089184 | 0.016251 | 0.154799 | 8.519661 | 0.577402 |
| 1 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.0085 | 0.001 | 13182.681372 | 842.416269 | 36.994473 | 45.335967 | 191.183969 | 5.621974 | 35.468645 | 6.189135 | 25.626771 | 4.542981 | 0.980962 | 0.262338 | 0.01262 | 0.11341 | 5.90972 | 0.568269 |
| 2 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 3 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.014 | 0.001 | 14347.472272 | 865.1912 | 38.924618 | 42.498292 | 245.347282 | 6.065827 | 36.549068 | 6.342976 | 29.559022 | 6.078583 | 1.004799 | 0.026049 | 0.012362 | 0.027869 | 6.134818 | 0.44031 |
| 4 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.006 | 0.001 | 11359.826166 | 398.808354 | 38.351538 | 35.11283 | 211.319577 | 5.625269 | 31.457117 | 5.143459 | 23.527645 | 12.577909 | 0.867042 | 0.039512 | 0.010175 | 0.016312 | 7.136776 | 0.511512 |
| 5 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.006 | 0.001 | 12756.337883 | 619.394774 | 50.671753 | 40.737497 | 375.535486 | 7.412915 | 41.57042 | 6.000838 | 27.053254 | 27.206014 | 1.024441 | 0.03496 | 0.012906 | 0.017003 | 11.635085 | 0.569886 |
| 6 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.006 | 0.001 | 15457.058054 | 899.500095 | 40.896228 | 38.923171 | 279.821146 | 7.684669 | 44.240518 | 8.256869 | 32.816672 | 11.471663 | 0.97767 | 0.077509 | 0.014367 | 0.018943 | 7.481523 | 0.573935 |
| 7 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.006 | 0.001 | 17890.662309 | 998.201546 | 45.404394 | 43.747367 | 427.318202 | 8.376312 | 46.478442 | 9.123958 | 35.343672 | 14.571893 | 1.142411 | 0.042442 | 0.017004 | 0.020126 | 8.038082 | 0.68719 |
| 8 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0075 | 0.001 | 17473.290669 | 905.622552 | 45.969994 | 49.265992 | 272.912776 | 7.716316 | 48.967096 | 8.753198 | 38.218362 | 7.303673 | 1.159622 | 0.069248 | 0.014839 | 0.021722 | 8.128869 | 0.700793 |
| 9 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.006 | 0.001 | 25411.662513 | 1709.858197 | 57.596293 | 63.408278 | 386.369163 | 10.613864 | 69.300728 | 13.314035 | 50.822242 | 9.160778 | 1.438946 | 0.037045 | 0.020781 | 0.028689 | 11.237761 | 0.705635 |
| 10 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 11 | Titiko - mud snails | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 12 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.1119 | 0.0048 | 9278.516506 | 10438.545347 | 26.122999 | 27.650273 | 191.842282 | 6.302344 | 37.37869 | 75.247891 | 49.645651 | 13.887514 | 3.036958 | 0.176368 | 0.052238 | 0.289213 | 5.886132 | 0.307983 |
| 13 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.1739 | 0.0083 | 3045.746777 | 9779.968577 | 9.266928 | 10.447118 | 87.293664 | 2.184617 | 15.890481 | 55.944349 | 32.241964 | 13.957774 | 2.454535 | 0.175301 | 0.040901 | 0.37774 | 1.59775 | 0.163822 |
| 14 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | 0.2742 | 0.031 | 1688.387959 | 8956.741669 | 7.28127 | 13.902505 | 168.717108 | 1.792196 | 8.471845 | 47.170575 | 27.372176 | 14.702333 | 2.223494 | 0.08714 | 0.02785 | 0.181421 | 1.095187 | 0.125211 |
| 15 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.4465 | 0.023 | 5834.559801 | 6688.965169 | 16.851306 | 17.802685 | 176.232594 | 4.467314 | 30.757799 | 37.809536 | 32.459531 | 13.655855 | 1.834363 | 0.20797 | 0.076393 | 0.216377 | 3.078978 | 0.251346 |
| 16 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.0452 | 0.0026 | 776.544571 | 4906.246349 | 2.118771 | 2.044623 | 94.758462 | 0.811912 | 3.536308 | 19.422672 | 10.663063 | 6.717349 | 1.087014 | 0.148339 | 0.017945 | 0.107886 | 0.532221 | 0.061831 |
| 17 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0309 | 0.0015 | 1498.423285 | 7484.46063 | 4.984812 | 7.125436 | 114.022396 | 1.854082 | 7.87016 | 50.283271 | 25.967419 | 14.104914 | 2.469336 | 0.28798 | 0.040168 | 0.312954 | 1.141666 | 0.112747 |
| 18 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.0251 | 0.0015 | 2461.750292 | 13737.544243 | 6.892372 | 7.129266 | 57.248232 | 1.750822 | 9.710237 | 62.373216 | 37.50133 | 25.564629 | 2.140799 | 0.451634 | 0.192853 | 0.09457 | 1.821799 | 0.203198 |
| 19 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.0276 | 0.0005 | 66.04758 | 12226.537642 | 0.39743 | 0.301674 | 15.563173 | 0.744519 | 5.384783 | 64.524262 | 29.971684 | 10.063169 | 2.498558 | 0.550304 | 0.033397 | 0.163524 | 0.352246 | 0.046921 |
| 20 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0774 | 0.0015 | 2677.861843 | 7411.325007 | 6.881487 | 7.085678 | 63.300903 | 1.397971 | 8.251588 | 37.786295 | 23.023848 | 10.121427 | 1.267161 | 0.047833 | 0.026337 | 0.141486 | 1.747061 | 0.108743 |
| 21 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.0431 | 0.0005 | 4030.26235 | 11537.183933 | 9.555477 | 9.870691 | 84.961609 | 1.806037 | 11.641405 | 48.033265 | 26.564529 | 9.444581 | 2.063936 | 0.179229 | 0.025591 | 0.21965 | 1.855676 | 0.136361 |
| 22 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 23 | Tuangi - cockles | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 24 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.0005 | 0.0005 | 787.639259 | 18087.464299 | 3.173642 | 3.555638 | 77.752051 | 2.079862 | 19.787828 | 6.039858 | 60.809418 | 30.604716 | 4.32808 | 0.091581 | 0.163688 | 0.088335 | 0.413524 | 0.2013 |
| 25 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.0005 | 0.0005 | 710.746254 | 17181.648407 | 2.634034 | 4.389376 | 20.956301 | 1.691787 | 15.865414 | 4.544055 | 49.425053 | 41.10754 | 3.770683 | 0.078458 | 0.155236 | 0.102417 | 0.437715 | 0.297476 |
| 26 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0005 | 0.0005 | 1166.470055 | 19532.038518 | 3.126026 | 3.65192 | 31.545677 | 1.952262 | 20.830599 | 5.549428 | 64.669416 | 23.676406 | 3.882477 | 0.184304 | 0.146906 | 0.077301 | 0.588409 | 0.151768 |
| 27 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 28 | Patiki - Flounder Liver | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 29 | Site 11 | 41˚ 15.680'S, 173˚ 5.875'E | 0.0317 | 0.0014 | 42.706757 | 12777.79052 | 0.46972 | 0.184253 | 10.980888 | 1.146335 | 0.522057 | 26.782493 | 136.388133 | 24.703859 | 5.558432 | 0.183492 | 0.377747 | 0.376774 | 0.073363 | 0.012373 |
| 30 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 31 | Patiki - Flounder Flesh | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 32 | Site 11 | 41˚ 15.680'S, 173˚ 5.875'E | 0.0042 | 0.0005 | 47.969064 | 10468.807791 | 0.173358 | 0.225523 | 1.891603 | 0.045363 | 0.225564 | 0.575889 | 37.250104 | 8.888757 | 1.288579 | 0.036508 | 0.00185 | 0.235958 | 0.035522 | 0.001545 |
| 33 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 34 | Patiki - Flounder Flesh | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 35 | NaN | NaN | NaN | NaN | 832.049226 | 9467.532397 | 2.155475 | 2.051407 | 48.205248 | 0.416026 | 2.524552 | 1.61799 | 81.990261 | 2.942643 | 1.165397 | 0.044176 | 0.007327 | 0.141226 | 0.308808 | 0.033755 |
| 36 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 37 | Kutai - green-lipped mussel | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | Nickel | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium |
| 38 | Site 12 | 41˚ 15.704'S, 173˚ 5.959'E | 0.003 | 0.0005 | 2382.523423 | 15487.941993 | 6.305763 | 6.132179 | 44.890628 | 1.771959 | 8.265882 | 5.052864 | 63.48187 | 13.710588 | 3.098001 | 0.040551 | 0.13091 | 0.157898 | 1.244749 | 0.14311 |
def parse_latlong(latlong):
if type(latlong) is str:
lat, long = re.findall(r"(\d+)˚ ([\d.]+)'(.)", latlong)
lat_deg, lat_min, lat_dir = lat
lon_deg, lon_min, lon_dir = long
latitude = (float(lat_deg) + float(lat_min) / 60) * (-1 if lat_dir == 'S' else 1)
longitude = (float(lon_deg) + float(lon_min) / 60) * (-1 if lon_dir == 'W' else 1)
return Point(longitude, latitude)
str(parse_latlong("41˚ 15.624'S, 173˚ 5.892'E"))
'POINT (173.0982 -41.2604)'
df["geometry"] = df.Location.apply(parse_latlong)
df = gpd.GeoDataFrame(df)
df.crs = 4326
df
| Sediment | Location | Total DDT Isomers | Dieldrin | Aluminuim | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | geometry | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.0075 | 0.001 | 16518.912655 | 1052.643599 | 42.463242 | 49.985737 | 273.474991 | 8.104016 | ... | 8.868234 | 35.506181 | 5.43789 | 0.991537 | 0.089184 | 0.016251 | 0.154799 | 8.519661 | 0.577402 | POINT (173.09820 -41.26040) |
| 1 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.0085 | 0.001 | 13182.681372 | 842.416269 | 36.994473 | 45.335967 | 191.183969 | 5.621974 | ... | 6.189135 | 25.626771 | 4.542981 | 0.980962 | 0.262338 | 0.01262 | 0.11341 | 5.90972 | 0.568269 | POINT (173.09822 -41.25953) |
| 2 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | POINT (173.09793 -41.25822) |
| 3 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.014 | 0.001 | 14347.472272 | 865.1912 | 38.924618 | 42.498292 | 245.347282 | 6.065827 | ... | 6.342976 | 29.559022 | 6.078583 | 1.004799 | 0.026049 | 0.012362 | 0.027869 | 6.134818 | 0.44031 | POINT (173.09652 -41.25838) |
| 4 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.006 | 0.001 | 11359.826166 | 398.808354 | 38.351538 | 35.11283 | 211.319577 | 5.625269 | ... | 5.143459 | 23.527645 | 12.577909 | 0.867042 | 0.039512 | 0.010175 | 0.016312 | 7.136776 | 0.511512 | POINT (173.09455 -41.25857) |
| 5 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.006 | 0.001 | 12756.337883 | 619.394774 | 50.671753 | 40.737497 | 375.535486 | 7.412915 | ... | 6.000838 | 27.053254 | 27.206014 | 1.024441 | 0.03496 | 0.012906 | 0.017003 | 11.635085 | 0.569886 | POINT (173.09288 -41.25860) |
| 6 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.006 | 0.001 | 15457.058054 | 899.500095 | 40.896228 | 38.923171 | 279.821146 | 7.684669 | ... | 8.256869 | 32.816672 | 11.471663 | 0.97767 | 0.077509 | 0.014367 | 0.018943 | 7.481523 | 0.573935 | POINT (173.08848 -41.25883) |
| 7 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.006 | 0.001 | 17890.662309 | 998.201546 | 45.404394 | 43.747367 | 427.318202 | 8.376312 | ... | 9.123958 | 35.343672 | 14.571893 | 1.142411 | 0.042442 | 0.017004 | 0.020126 | 8.038082 | 0.68719 | POINT (173.08700 -41.25658) |
| 8 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0075 | 0.001 | 17473.290669 | 905.622552 | 45.969994 | 49.265992 | 272.912776 | 7.716316 | ... | 8.753198 | 38.218362 | 7.303673 | 1.159622 | 0.069248 | 0.014839 | 0.021722 | 8.128869 | 0.700793 | POINT (173.09685 -41.25855) |
| 9 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.006 | 0.001 | 25411.662513 | 1709.858197 | 57.596293 | 63.408278 | 386.369163 | 10.613864 | ... | 13.314035 | 50.822242 | 9.160778 | 1.438946 | 0.037045 | 0.020781 | 0.028689 | 11.237761 | 0.705635 | POINT (173.09020 -41.25885) |
| 10 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 11 | Titiko - mud snails | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 12 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.1119 | 0.0048 | 9278.516506 | 10438.545347 | 26.122999 | 27.650273 | 191.842282 | 6.302344 | ... | 75.247891 | 49.645651 | 13.887514 | 3.036958 | 0.176368 | 0.052238 | 0.289213 | 5.886132 | 0.307983 | POINT (173.09820 -41.26040) |
| 13 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.1739 | 0.0083 | 3045.746777 | 9779.968577 | 9.266928 | 10.447118 | 87.293664 | 2.184617 | ... | 55.944349 | 32.241964 | 13.957774 | 2.454535 | 0.175301 | 0.040901 | 0.37774 | 1.59775 | 0.163822 | POINT (173.09822 -41.25953) |
| 14 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | 0.2742 | 0.031 | 1688.387959 | 8956.741669 | 7.28127 | 13.902505 | 168.717108 | 1.792196 | ... | 47.170575 | 27.372176 | 14.702333 | 2.223494 | 0.08714 | 0.02785 | 0.181421 | 1.095187 | 0.125211 | POINT (173.09793 -41.25822) |
| 15 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.4465 | 0.023 | 5834.559801 | 6688.965169 | 16.851306 | 17.802685 | 176.232594 | 4.467314 | ... | 37.809536 | 32.459531 | 13.655855 | 1.834363 | 0.20797 | 0.076393 | 0.216377 | 3.078978 | 0.251346 | POINT (173.09652 -41.25838) |
| 16 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.0452 | 0.0026 | 776.544571 | 4906.246349 | 2.118771 | 2.044623 | 94.758462 | 0.811912 | ... | 19.422672 | 10.663063 | 6.717349 | 1.087014 | 0.148339 | 0.017945 | 0.107886 | 0.532221 | 0.061831 | POINT (173.09455 -41.25857) |
| 17 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0309 | 0.0015 | 1498.423285 | 7484.46063 | 4.984812 | 7.125436 | 114.022396 | 1.854082 | ... | 50.283271 | 25.967419 | 14.104914 | 2.469336 | 0.28798 | 0.040168 | 0.312954 | 1.141666 | 0.112747 | POINT (173.09288 -41.25860) |
| 18 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.0251 | 0.0015 | 2461.750292 | 13737.544243 | 6.892372 | 7.129266 | 57.248232 | 1.750822 | ... | 62.373216 | 37.50133 | 25.564629 | 2.140799 | 0.451634 | 0.192853 | 0.09457 | 1.821799 | 0.203198 | POINT (173.08848 -41.25883) |
| 19 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.0276 | 0.0005 | 66.04758 | 12226.537642 | 0.39743 | 0.301674 | 15.563173 | 0.744519 | ... | 64.524262 | 29.971684 | 10.063169 | 2.498558 | 0.550304 | 0.033397 | 0.163524 | 0.352246 | 0.046921 | POINT (173.08700 -41.25658) |
| 20 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0774 | 0.0015 | 2677.861843 | 7411.325007 | 6.881487 | 7.085678 | 63.300903 | 1.397971 | ... | 37.786295 | 23.023848 | 10.121427 | 1.267161 | 0.047833 | 0.026337 | 0.141486 | 1.747061 | 0.108743 | POINT (173.09685 -41.25855) |
| 21 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.0431 | 0.0005 | 4030.26235 | 11537.183933 | 9.555477 | 9.870691 | 84.961609 | 1.806037 | ... | 48.033265 | 26.564529 | 9.444581 | 2.063936 | 0.179229 | 0.025591 | 0.21965 | 1.855676 | 0.136361 | POINT (173.09020 -41.25885) |
| 22 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 23 | Tuangi - cockles | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 24 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.0005 | 0.0005 | 787.639259 | 18087.464299 | 3.173642 | 3.555638 | 77.752051 | 2.079862 | ... | 6.039858 | 60.809418 | 30.604716 | 4.32808 | 0.091581 | 0.163688 | 0.088335 | 0.413524 | 0.2013 | POINT (173.09820 -41.26040) |
| 25 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.0005 | 0.0005 | 710.746254 | 17181.648407 | 2.634034 | 4.389376 | 20.956301 | 1.691787 | ... | 4.544055 | 49.425053 | 41.10754 | 3.770683 | 0.078458 | 0.155236 | 0.102417 | 0.437715 | 0.297476 | POINT (173.09822 -41.25953) |
| 26 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0005 | 0.0005 | 1166.470055 | 19532.038518 | 3.126026 | 3.65192 | 31.545677 | 1.952262 | ... | 5.549428 | 64.669416 | 23.676406 | 3.882477 | 0.184304 | 0.146906 | 0.077301 | 0.588409 | 0.151768 | POINT (173.09288 -41.25860) |
| 27 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 28 | Patiki - Flounder Liver | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 29 | Site 11 | 41˚ 15.680'S, 173˚ 5.875'E | 0.0317 | 0.0014 | 42.706757 | 12777.79052 | 0.46972 | 0.184253 | 10.980888 | 1.146335 | ... | 26.782493 | 136.388133 | 24.703859 | 5.558432 | 0.183492 | 0.377747 | 0.376774 | 0.073363 | 0.012373 | POINT (173.09792 -41.26133) |
| 30 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 31 | Patiki - Flounder Flesh | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 32 | Site 11 | 41˚ 15.680'S, 173˚ 5.875'E | 0.0042 | 0.0005 | 47.969064 | 10468.807791 | 0.173358 | 0.225523 | 1.891603 | 0.045363 | ... | 0.575889 | 37.250104 | 8.888757 | 1.288579 | 0.036508 | 0.00185 | 0.235958 | 0.035522 | 0.001545 | POINT (173.09792 -41.26133) |
| 33 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 34 | Patiki - Flounder Flesh | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 35 | NaN | NaN | NaN | NaN | 832.049226 | 9467.532397 | 2.155475 | 2.051407 | 48.205248 | 0.416026 | ... | 1.61799 | 81.990261 | 2.942643 | 1.165397 | 0.044176 | 0.007327 | 0.141226 | 0.308808 | 0.033755 | None |
| 36 | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | None |
| 37 | Kutai - green-lipped mussel | NaN | Total DDT Isomers | Dieldrin | Aluminium | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | None |
| 38 | Site 12 | 41˚ 15.704'S, 173˚ 5.959'E | 0.003 | 0.0005 | 2382.523423 | 15487.941993 | 6.305763 | 6.132179 | 44.890628 | 1.771959 | ... | 5.052864 | 63.48187 | 13.710588 | 3.098001 | 0.040551 | 0.13091 | 0.157898 | 1.244749 | 0.14311 | POINT (173.09932 -41.26173) |
39 rows × 21 columns
float_cols = df.columns[2:-1]
float_cols
Index(['Total DDT Isomers', 'Dieldrin', 'Aluminuim', 'Sulphur', 'Vanadium',
'Chromium', 'Manganese', 'Cobalt', 'Nickel', 'Copper', 'Zinc',
'Arsenic', 'Selenium', 'Silver', 'Cadmium', 'Mercury', 'Lead',
'Uranium'],
dtype='object')
sediment = df.iloc[:10].copy()
sediment[float_cols] = sediment[float_cols].astype(float)
sediment
| Sediment | Location | Total DDT Isomers | Dieldrin | Aluminuim | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | geometry | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.0075 | 0.001 | 16518.912655 | 1052.643599 | 42.463242 | 49.985737 | 273.474991 | 8.104016 | ... | 8.868234 | 35.506181 | 5.437890 | 0.991537 | 0.089184 | 0.016251 | 0.154799 | 8.519661 | 0.577402 | POINT (173.09820 -41.26040) |
| 1 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.0085 | 0.001 | 13182.681372 | 842.416269 | 36.994473 | 45.335967 | 191.183969 | 5.621974 | ... | 6.189135 | 25.626771 | 4.542981 | 0.980962 | 0.262338 | 0.012620 | 0.113410 | 5.909720 | 0.568269 | POINT (173.09822 -41.25953) |
| 2 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | POINT (173.09793 -41.25822) |
| 3 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.0140 | 0.001 | 14347.472272 | 865.191200 | 38.924618 | 42.498292 | 245.347282 | 6.065827 | ... | 6.342976 | 29.559022 | 6.078583 | 1.004799 | 0.026049 | 0.012362 | 0.027869 | 6.134818 | 0.440310 | POINT (173.09652 -41.25838) |
| 4 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.0060 | 0.001 | 11359.826166 | 398.808354 | 38.351538 | 35.112830 | 211.319577 | 5.625269 | ... | 5.143459 | 23.527645 | 12.577909 | 0.867042 | 0.039512 | 0.010175 | 0.016312 | 7.136776 | 0.511512 | POINT (173.09455 -41.25857) |
| 5 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0060 | 0.001 | 12756.337883 | 619.394774 | 50.671753 | 40.737497 | 375.535486 | 7.412915 | ... | 6.000838 | 27.053254 | 27.206014 | 1.024441 | 0.034960 | 0.012906 | 0.017003 | 11.635085 | 0.569886 | POINT (173.09288 -41.25860) |
| 6 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.0060 | 0.001 | 15457.058054 | 899.500095 | 40.896228 | 38.923171 | 279.821146 | 7.684669 | ... | 8.256869 | 32.816672 | 11.471663 | 0.977670 | 0.077509 | 0.014367 | 0.018943 | 7.481523 | 0.573935 | POINT (173.08848 -41.25883) |
| 7 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.0060 | 0.001 | 17890.662309 | 998.201546 | 45.404394 | 43.747367 | 427.318202 | 8.376312 | ... | 9.123958 | 35.343672 | 14.571893 | 1.142411 | 0.042442 | 0.017004 | 0.020126 | 8.038082 | 0.687190 | POINT (173.08700 -41.25658) |
| 8 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0075 | 0.001 | 17473.290669 | 905.622552 | 45.969994 | 49.265992 | 272.912776 | 7.716316 | ... | 8.753198 | 38.218362 | 7.303673 | 1.159622 | 0.069248 | 0.014839 | 0.021722 | 8.128869 | 0.700793 | POINT (173.09685 -41.25855) |
| 9 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.0060 | 0.001 | 25411.662513 | 1709.858197 | 57.596293 | 63.408278 | 386.369163 | 10.613864 | ... | 13.314035 | 50.822242 | 9.160778 | 1.438946 | 0.037045 | 0.020781 | 0.028689 | 11.237761 | 0.705635 | POINT (173.09020 -41.25885) |
10 rows × 21 columns
for col in tqdm(float_cols):
ax = sediment.plot(col, figsize=(10,10), legend=True, legend_kwds={"orientation": "horizontal", "pad": 0.05}, cmap="Reds", markersize = 100, edgecolor = "black")
ax.set_title(f"Sediment {col} in Mapua")
ax.ticklabel_format(useOffset=False, style='plain')
xmin, ymin, xmax, ymax = sediment.total_bounds
pad = .002
ax.set_xlim(xmin-pad, xmax+pad)
ax.set_ylim(ymin-pad, ymax+pad)
cx.add_basemap(ax, source=cx.providers.Esri.WorldImagery, crs=sediment.crs, zoom=18)
file_name = f"Sediment {col} in Mapua.png"
ax.figure.savefig(file_name, bbox_inches="tight")
0%| | 0/18 [00:00<?, ?it/s]
titiko = df.iloc[12:22].copy()
titiko[float_cols] = titiko[float_cols].astype(float)
titiko
| Sediment | Location | Total DDT Isomers | Dieldrin | Aluminuim | Sulphur | Vanadium | Chromium | Manganese | Cobalt | ... | Copper | Zinc | Arsenic | Selenium | Silver | Cadmium | Mercury | Lead | Uranium | geometry | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 12 | Site 1 | 41˚ 15.624'S, 173˚ 5.892'E | 0.1119 | 0.0048 | 9278.516506 | 10438.545347 | 26.122999 | 27.650273 | 191.842282 | 6.302344 | ... | 75.247891 | 49.645651 | 13.887514 | 3.036958 | 0.176368 | 0.052238 | 0.289213 | 5.886132 | 0.307983 | POINT (173.09820 -41.26040) |
| 13 | Site 2 | 41˚ 15.572'S, 173˚ 5.893'E | 0.1739 | 0.0083 | 3045.746777 | 9779.968577 | 9.266928 | 10.447118 | 87.293664 | 2.184617 | ... | 55.944349 | 32.241964 | 13.957774 | 2.454535 | 0.175301 | 0.040901 | 0.377740 | 1.597750 | 0.163822 | POINT (173.09822 -41.25953) |
| 14 | Site 3 | 41˚ 15.493'S, 173˚ 5.876'E | 0.2742 | 0.0310 | 1688.387959 | 8956.741669 | 7.281270 | 13.902505 | 168.717108 | 1.792196 | ... | 47.170575 | 27.372176 | 14.702333 | 2.223494 | 0.087140 | 0.027850 | 0.181421 | 1.095187 | 0.125211 | POINT (173.09793 -41.25822) |
| 15 | Site 4 | 41˚ 15.503'S, 173˚ 5.791'E | 0.4465 | 0.0230 | 5834.559801 | 6688.965169 | 16.851306 | 17.802685 | 176.232594 | 4.467314 | ... | 37.809536 | 32.459531 | 13.655855 | 1.834363 | 0.207970 | 0.076393 | 0.216377 | 3.078978 | 0.251346 | POINT (173.09652 -41.25838) |
| 16 | Site 5 | 41˚ 15.514'S, 173˚ 5.673'E | 0.0452 | 0.0026 | 776.544571 | 4906.246349 | 2.118771 | 2.044623 | 94.758462 | 0.811912 | ... | 19.422672 | 10.663063 | 6.717349 | 1.087014 | 0.148339 | 0.017945 | 0.107886 | 0.532221 | 0.061831 | POINT (173.09455 -41.25857) |
| 17 | Site 6 | 41˚ 15.516'S, 173˚ 5.573'E | 0.0309 | 0.0015 | 1498.423285 | 7484.460630 | 4.984812 | 7.125436 | 114.022396 | 1.854082 | ... | 50.283271 | 25.967419 | 14.104914 | 2.469336 | 0.287980 | 0.040168 | 0.312954 | 1.141666 | 0.112747 | POINT (173.09288 -41.25860) |
| 18 | Site 7 | 41˚ 15.530'S, 173˚ 5.309'E | 0.0251 | 0.0015 | 2461.750292 | 13737.544243 | 6.892372 | 7.129266 | 57.248232 | 1.750822 | ... | 62.373216 | 37.501330 | 25.564629 | 2.140799 | 0.451634 | 0.192853 | 0.094570 | 1.821799 | 0.203198 | POINT (173.08848 -41.25883) |
| 19 | Site 8 | 41˚ 15.395'S, 173˚ 5.220'E | 0.0276 | 0.0005 | 66.047580 | 12226.537642 | 0.397430 | 0.301674 | 15.563173 | 0.744519 | ... | 64.524262 | 29.971684 | 10.063169 | 2.498558 | 0.550304 | 0.033397 | 0.163524 | 0.352246 | 0.046921 | POINT (173.08700 -41.25658) |
| 20 | Site 9 | 41˚ 15.513'S, 173˚ 5.811'E | 0.0774 | 0.0015 | 2677.861843 | 7411.325007 | 6.881487 | 7.085678 | 63.300903 | 1.397971 | ... | 37.786295 | 23.023848 | 10.121427 | 1.267161 | 0.047833 | 0.026337 | 0.141486 | 1.747061 | 0.108743 | POINT (173.09685 -41.25855) |
| 21 | Site 10 | 41˚ 15.531'S, 173˚ 5.412'E | 0.0431 | 0.0005 | 4030.262350 | 11537.183933 | 9.555477 | 9.870691 | 84.961609 | 1.806037 | ... | 48.033265 | 26.564529 | 9.444581 | 2.063936 | 0.179229 | 0.025591 | 0.219650 | 1.855676 | 0.136361 | POINT (173.09020 -41.25885) |
10 rows × 21 columns
for col in tqdm(float_cols):
ax = titiko.plot(col, figsize=(10,10), legend=True, legend_kwds={"orientation": "horizontal", "pad": 0.05}, cmap="Reds", markersize = 100, edgecolor = "black")
ax.set_title(f"Tītiko (mud-flat snails) {col} in Mapua")
ax.ticklabel_format(useOffset=False, style='plain')
xmin, ymin, xmax, ymax = titiko.total_bounds
pad = .002
ax.set_xlim(xmin-pad, xmax+pad)
ax.set_ylim(ymin-pad, ymax+pad)
cx.add_basemap(ax, source=cx.providers.Esri.WorldImagery, crs=sediment.crs, zoom=18)
file_name = f"Tītiko (mud-flat snails) {col} in Mapua.png"
ax.figure.savefig(file_name, bbox_inches="tight")
0%| | 0/18 [00:00<?, ?it/s]